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Abstract. We discuss topological aspects of electronic properties of graphene, in- 
cluding edge effects, with the tight-binding model on a honeycomb lattice and 
its extensions to show the following: (i) Appearance of the pair of massless Dirac 
dispersions, which is the origin of anomalous properties including a peculiar quan- 
tum Hall effect (QHE), is not accidental to honeycomb, but is rather generic for 
a class of two-dimensional lattices that interpolate between square and 7r-flux lat- 
tices. Persistence of the peculiar QHE is interpreted as a topological stability, (ii) 
While we have the massless Dirac dispersion only around E = 0, the anomalous 
QHE associated with the Dirac cone unexpectedly persists for a wide range of 
the chemical potential. The range is bounded by van Hove singularities, at which 
we predict a transition to the ordinary fermion behavior accompanied by huge 
jumps in the QHE with a sign change, (iii) For edges we establish a coincidence 
between the quantum Hall effect in the bulk and the quantum Hall effect for the 
edge states, which is a manifestation of the topological bulk-edge correspondence. 
We have also explicitly shown that the E — edge states in honeycomb in zero 
magnetic field persist in magnetic field. 

1 Introduction 

One fascinating aspect of the condensed-matter physics is that we can have various field theories 
effectively realized on low-energy scales. Electronic structure in a honeycomb lattice is of partic- 
ular interest, since the system has zero-mass Dirac cone dispersion at Brillouin zone corners, so 
we have got an intriguing mixture of a continuum limit and the lattice effect. Graphene, a single 
layer carbon atoms, provides a prototypical realization of characteristic electronic structures of 
a honeycomb lattice. This provides interesting problems in condensed matter physics especially 
in its topological aspects, which is the subject-matter of the present article. Its band dispersion 
is composed of a pair of fc-linear bands that touch with each other around E = 0, so that the 
graphene is a condensed-matter realization of zero-mass Dirac fermions around E = 0. Now, 
the Dirac particles have been known to provide topologically non-trivial structures |l|2j . Physics 
on the honeycomb lattice thus provides potential links between the field-theoretic concepts and 
various condensed matter phenomena [31415) . 

Indeed, recent experiments have revealed a peculiar type of quantum Hall effect [BT7I8) . 
Graphite itself has been studied extensively from many different points of view. There has 
also been a substantial amount of development as basic materials in nano-technologies [91 1 0) . 
where a special interest is characteristic boundary magnetic properties which originate from 
edge states [Tl|12| . The fact that the honeycomb lattice as interesting edge properties as well 
as interesting bulk properties is no accident, but this is in fact a realization of hidden physical 
structures as a topological order |13|14H5j . 
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Topologically non-trivial ground states are, unlike ordinary symmetry-broken states, char- 
acterized by geometrical phases rather than by local order parameters. Geometrical phases 
which are fundamental in quantum physics provide basic tools for characterizing topological 
orders [T4I15I16] . Appearance of the edge states is one of the important physical consequences 
of the non-trivial topological orders [TTTl 8 19 20 21j. Physical properties of the bulk go hand 
in hand with those of the edge in the topologically ordered states. So topological properties, 
which may not so apparent in the bulk, may emerge around the edges or defects with physical 
consequences. |17|18|19l20] . This bulk-edge correspondence is shared by various topologically 
ordered states, where a typical one is the quantized Hall (QH) effects [19 20 . Among various 
topologically ordered states, edge states of the Dirac particle is interesting, especially when the 
particle is massless l22l23| . This is in fact directly related with the edge states in graphene. To 
be more precise, honeycomb has a bipartite lattice structure (i.e., a chiral symmetry), which 
has a fundamental importance in the existence of the E = edge states. From the topological 
point of view, appearance of the magnetic moments in the edge states in graphene is understood 
topologically as a chiral symmetry breaking in edge states }24I25| . 

In this article, we discuss results obtained in the reference [26] focusing on the topological 
aspects of graphene. 



2 Topological Stability of Dirac Fermions 

We consider a tight-binding Hamiltonian, H, for a honeycomb lattice with nearest-neighbor 
hopping H' 

H ^honeycomb H 
^honeycomb = t [ C I ( J ) Co (j ) + e* 2 *** cj ( j )c Q ( j - e 2 ) + c\{j + ei)c(i)] + h.C, 

3 

H'=t'J2 e- i2 ^ +1 ^cl(j + e x - e 2 )c (j) + h.c. 
3 

where t = — 1 is the transfer energy (taken to be the unit of energy). Since the honeycomb, a 
non-Bravais lattice with two sites per unit cell is bipartite, we have denned two kinds of fermion 
operators, c (j) and c.(j) as in Fig. [1] where j = j%ei + labels the unit cell with two 
translation vectors e\ = (3/2, v3/2)a and e 2 = (0, \/?>)a. 
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Fig. 2. (a) Locus (each circle) of A(k) is shown on the complex Zi-plane when we vary k x : — > 2n 
for several fixed values of k y for t' /t — —1/2. The red circle traverses the origin, which implies the 
appearance of the gapless Dirac fermion. The doubling of the Dirac fermion is also clear from this 
geometric consideration, (b) Corresponding energy dispersion for t' /t = —1/2. 



We then apply an external magnetic field, which is characterized by 

ct> = BS 6 /(27r)=p/q, 

the magnetic flux penetrating each hexagon of area Sq = (3-\/3/2)a 2 . 

Let us first discuss the energy dispersion in zero magnetic field. As is well known, the 
spectrum of the Hamiltonian -ffhoneycomb consists of two bands touching with each other at 
two points (K and K') in the Brillouin zone, which are understood as doubled Dirac fermions. 
Around the zero-gap points the energy dispersion is generically linear in (k x , k y ) space. 

Although the appearance of the massless Dirac fermions may seem just an accident for 
the honeycomb lattice, we shall show here that the appearance of massless Dirac fermions 
is, surprisingly enough, a generic property and stable against e.g. introduction of a diagonal 
hopping H' if we preserve the chiral symmetry (bipartite lattice with every transfer between A 
and B sites). In momentum representation, the Hamiltonian in zero magnetic field is given as 

H = J ^cHk)h{k)c{k) 

«*> "(SSI)' 
h(k)-t( ° ^ (fc) ^l 

A(k) =1 + e"* 2 + e~ ifcl [1 + {t'/t)e ik2 ] . 
The chiral symmetry is expressed as 

{h,i} = 0, J = (T z= (l > 

which implies that the energy bands are symmetric about E — 0. For the honeycomb lattice 
A on the complex plane delineates a circle centered at Co = 1 + e lk2 with a radius r = 

yl + 2(t'/t) cos &2 + (t' /t) 2 when k\ is varied from to 2ir for a fixed value of k y (Fig. 2). By 
examining Co and r, we can show that we have a Dirac cone when the circle goes through the 
origin, which occurs for — 3 < t'/t < 1. This shows a topological stability for the appearance of 
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the Dirac fermions. That is, Dirac cone persists when the lattice deviates from the honeycomb, 
even up to the limit of the square lattice where t! jt = 1 in the present notation, unless the 
chiral symmetry is not broken. We can extend the present line of consideration to the general 
condition for the appearance of Dirac cones, which will be published elsewhere |27j. 

3 Hall conductance of the graphene as the Chern number 

The Hall conductivity of the noninteracting 2D electron systems, as described with the linear 
response theory. [28] may be regarded as a topological quantum number .[29 Namely, when the 
Fermi energy lies in the n-th gap, the Hall conductivity is given by 

e 2 . , 

0~xy = J^Cf[Ef), 

where cf is an integer called the Chern number, which describes how the wave function responds 
to a vector potential generated by Berry's gauge potential in the Brillouin zone. If we adopt a 
Bloch function in the j-th energy band <p 3 (fe) assuming that the energy gap opens over the entire 
Brillouin zone at each gap in the sequence ei(fe) < • • ■ < € n {k) < Ep < e„+i(fc) < • • • < e2 q (k). 
We can then define a non-Abelian Berry's gauge potential, 

A = 

Here A(k) is an n x n matrix-valued one-form, which is a difference from the ordinary, scalar 
Berry's gauge potential. Then the Chern number is given as the U(l) part of the U(n) gauge 
potential, 

cf(E f ) = J Tr dA. 

This formulation, a non-Abelian extension of Berry's gauge potential. [14|15j may seem too 
elaborate, but is useful when there are multiple bands below the Fermi energy. Namely, this 
formula holds even if some gaps in the Fermi sea are closed, as long as the bands in question, 
e n (k) and e n +i(k), do not overlap. In the special case when all the bands in the Fermi sea are 
separated with each other, the Chern number, Eq.j^Q), is simply the sum of the Chern numbers 
assigned to individual bands. 

To evaluate the expression, here we have employed a technique developed recently in the 
lattice gauge theory by two of the present authors [30|31j . The numerical results for the Chern 
number cp, i.e., the Hall conductivity per spin, as a function of the chemical potential Ep for 
honeycomb (t 1 = 0) is shown in Fig. 3. While the integral Chern number cp is defined only for 
each gap, we have plotted Cf as step functions, which makes sense as long as the magnetic field 
is not too large, as in Fig. [3]with <f> = 1/31, since the width of each Landau band is then much 
smaller than the size of gaps. 

The result shows a striking answer to a key question: what is the fate of the Dirac-like 
behavior as we go away from E = 0. When the Fermi energy is swept in the honeycomb lattice, 
the Dirac-like Hall conductivity with steps of two (or four when spin degeneracy is included) 
in units of —e 2 /h exist around E = as has been noted by many authors. |6|7|8|32j Let us 
call this Dirac-Landau behavior. Now, we can see that this Dirac-Landau behavior persists all 
the way up to some energy, E v g — \t\, which we can identify to be the van Hove singularities 
appearing in Fig. [3l We can notice that huge steps accompanied by a sign change in the Hall 
conductivity occurs at these points. This result implies the following: The effective theory near 
the zero energy is Dirac-like fermions in the continuum limit, for which an unconventional 
quantization of the Hall conductivity has been derived in the continuum limit from the Dirac 
Landau levels. [6l9] The present calculation leads to the conclusion that the unusual property 
extends to an unexpectedly wide region of energy in the actual lattice fermion system. 
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Fig. 3. (a) The Chern number (Hall conductivity per spin in units of —e 2 /h) for magnetic field 
<f> = 1/31 is plotted against the Fermi energy Ef for t' /t = 0. The red lines indicate the steps of two 
in the Chern number sequence (<r X y = ±(2iV + l)e 2 /h, N: integer), while the blue lines step of one 
{<Jxy = ±iVe 2 /h). (b) Density of states for the honeycomb lattice. 



Outside the van Hove energies (i.e., in the band-edge regions), we recover the conventional 
quantum Hall effect (QHE) where the step changes by one in units of —e 2 /h (which we call 
Fermi-Landau). This automatically implies that a huge step accompanied by a sign change 
has to occur at the boundary between Dirac-Landau and Fermi-Landau behaviors. In other 
words the bands situated at the van-Hove singularities induce a change, which is topological 
in that the relevant quantum numbers are topological. The Hall conductivity per spin on the 
honeycomb lattice for the entire energy region is summarized as ( N = 0, 1, 2, • • ■ ) 

{+(N + 1) E-p < —E v h, ■ Landau index from the band bottom 

— (27V + 1) —E v h < Ep < 0, : Dirac-Landau index from the zero energy 

+(2JV + 1) < Ep < +E v h, ■ Dirac-Landau index from the zero energy 

— (TV + 1) E v h < Ef, : Landau index from the band top 

We should again double these numbers if we include the spin degeneracy. While the unconven- 
tional quantization around the zero energy has been beautifully observed experimentally [7 8J, 
pushing the chemical potential further to approach the van Hove energies should detect the 
topological phase transition. 



4 Bulk-Edge Correspondence in Graphene: Edge states of the Dirac 
Fermions 

When a system has a nontrivial topological structure, edge states should generically exhibit 
characteristic properties. For the honeycomb lattice, edge states and their flat dispersion have 
been intensively discussed with its relevance to spin alignment [1 1 1 1 2 j . Another way of saying 
is that the honeycomb system has a bipartite symmetry, which guarantees the presence of the 
zero mode edge states[24], i.e., there exist macroscopic edge states at the zero energy unless 
the bipartite symmetry is broken. This gives rise to dispersionless edge states. The fiat bands 
are unstable against bipartite symmetry breaking (Peierls instability), which is, e.g., realized 
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Fig. 4. The energy spectra for the cylindrical system in Fig. [T] with the flux (j> = 1/5 for t'/t = 
(honeycomb). k y is a momentum along the cylindrical direction. Blue lines represent edge states 
localized at the right (bearded) edge, while red ones the left (zigzag) edge, with the bulk energy bands 
shown as grey regions. 



when an antiferromagnetic spin ordering occurs. This may be viewed as a topological origin of 
the boundary spin moments. [25] 

Now, the problem of edge states vs bulk states is particularly interesting for the QHE, 
since the problem of how the Dirac-Landau QHE oc (2iV + 1) is related to the edge states is 
fundamentally interesting. One of the present authors has shown, for 2D square (or anisotropic 
square) systems, that the edge states whose energy dispersions run across the Landau bands 
have topological QHE numbers, for which cr^ ge = cr~y by identifying the connection between 
the topological integers for the bulk and for the edge states. [20121 j So let us now focus on the 
edge states of the honeycomb lattice in terms of the topological structure. 

We follow the analysis in refs 19 20 21j to reduce the system to a one-dimensional model by 
a partial Fourier transform along the second direction (j'2) of the fermion operators, c ,»(J) = 

Co,»(ji, fe). This yields a ^-dependent series of one-dimensional Hamiltonian 



rl/2 v 



Jk 2 j2 



k 2 

Hm{k 2 ) = ^2 [*o.0'i,fe)40'i, fc 2)c.(ji,fc 2 ) +t.o(Ji,k 2 )cl(ji + 1, k 2 )c (ji, k 2 ) 



h.c. 



31 



where /^-dependent hopping parameters are 

to.(jiM) =t(l + e ifc2 - Wl ), 
t.o{jiM =t [l + (t'/t) e ^-^(n+i/2) 



We consider cylindrical systems with a zigzag edge at one end and a bearded (or Klein's) 
edge at the other. We have adopted this geometry, since zigzag-bearded system has an E = 
edge state in zero magnetic field that extend over the entire Brillouin zone, as first pointed out 
in l33l. 
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Fig. 5. The charge density of the edge states at the zigzag edge for <f> = 1/3. (a) £7 = —0.95 and 
E = -1.5. 



A remarkable feature in the spectrum shown here is that, on top of the edge states across 
adjacent Landau bands, there always exist exactly E = edge states with a flat dispersion, 
which has to exist since the number of the edge states are odd, 2q — 1 . Namely, a system with 
zigzag edges has a zero-energy edge mode in some region of ki , a system with bearded edges in 
the other region, and the present system has the mode over the entire region. This is a simple 
extension of the discussion in the absence of a magnetic field. [24l25j 

With Laughlin's argument [17] . supplemented by the behavior of the edge state (as a function 
of the momentum fo), one can identify a topological number I(Ep), the number of electrons 
carried by Laughlin's adiabatic procedure, for each edge branch when the chemical potential 
E-p traverses the branch. |17|19|20|21] As shown in Fig. 4, we can readily confirmed a bulk-edge 
correspondence, 

c F (£ F ) = I{E F ) 

bulk _ edge 
xy xy 

for Dirac fermions in graphene. 

As for the edge states, we display in Fig. 5 typical wave functions in real space. We can see 
that each of the left and right edge states decays into the bulk, where the decay length depends 
on whether its energy is close to the touching of the edge mode and the Landau level. 
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5 summary 

In summary, we have shown that the Landau levels are divided by the van Hove singularities 
into two regimes: One is effectively described by Dirac particles, and the other by ordinary 
finite-mass fermions. We have confirmed this both from the Chern number argument and by 
the bulk-edge correspondence. Topological considerations are useful throughout. 

The properties revealed here for the whole spectrum of the honeycomb lattice may be 
experimentally observable if the chemical potential can be varied over a wide range. In real 
graphene samples there may be disorder, in which case we are talking about a dirty Hofstadter 
system. However, non-monotonic behaviors should survive the disorder as far as the degree of 
disorder is not too large, as has been indicated by a numerical calculation for a dirty Hofstadter 
system [34). 
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